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Abstract 

Background: Tine uniform delivery of laser energy is particularly important for safe 
and effective laser hair removal (LHR) treatment. Although it is necessary to 
quantitatively assess the spatial distribution of the delivered laser, laser spots are 
difficult to trace owing to a lack of visual cues. This study proposes a novel preclinic 
tool to evaluate operator proficiency in LHR treatment and applies this tool to train 
novice operators and compare two different treatment techniques (sliding versus 
spot- by-spot). 

Methods: A simulation bed is constructed to visualize the irradiated laser spots. Six 
novice operators are recruited to perform four sessions of simulation while changing 
the treatment techniques and the presence of feedback (sliding without feedback, 
sliding with feedback, spot-by-spot without feedback, and spot-by-spot with feedback). 
Laser distribution maps (LDMs) are reconstructed through a series of images processed 
from the recorded video for each simulation session. Then, an experienced dermatologist 
classifies the collected LDMs into three different performance groups, which are 
quantitatively analyzed in terms of four performance indices. 

Results: The performance groups are characterized by using a combination of four 
proposed indices. The best-performing group exhibited the lowest amount of randomness 
in laser delivery and accurate estimation of mean spot distances. The training was 
only effective in the sliding treatment technique. After the training, omission errors 
decreased by 6.32% and better estimation of the mean spot distance of the actual size 
of the laser-emitting window was achieved. Gels required operators to be trained when 
the spot-by-spot technique was used, and imposed difficulties in maintaining regular 
laser delivery when the sliding technique was used. 

Conclusions: Because the proposed system is simple and highly affordable, it is 
expected to benefit many operators in clinics to train and maintain skilled performance 
in LHR treatment, which will eventually lead to accomplishing a uniform laser delivery 
for safe and effective LHR treatment. 

Keywords: Laser visualization. Performance evaluation, Photomedicine, Simulation bed. 
Training tool 
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Background 

Approximately two decades after the emergence of the first FDA-approved laser, laser 
hair removal (LHR) treatment has become one of the most successful applications of 
lasers in medicine [1]. According to recent statistics from the American Society for 
Aesthetic Plastic Surgery (AS APS), more than 1.2 million LHR procedures were per- 
formed in the U.S. during 2012, which were the third most frequent treatments in all 
cosmetic procedures and the second most frequent treatments for men [2]. The preva- 
lence of LHR is assumed to be largely attributed to increased societal concerns about 
aesthetics, as well as its proven safety, efficacy, and greater simplicity than conventional 
epilation methods [3-5]. 

The idea of LHR is based on the theory of selective photothermolysis. This elaborates 
the differences in absorption rates of light energy between the hair follicle and tissue 
owing to the differences in chromophores [6]. Therefore, hair follicles can be selectively 
destroyed by exposing the laser to the target area without aiming at each follicle. Ac- 
cording to more recent findings, the essence of photoepilation is not just the heat gen- 
eration in the hair follicle but the conduction of heat to the hair stem cells, because the 
hair stem cell itself lacks an appreciable amount of chromophores and is located out- 
side the outer root sheath of the hair follicle [7]. Therefore, the delivery of an appropri- 
ate dose of laser light is paramount for safe and effective photoepilation. 

Failure in choosing the correct laser intensity level often causes side effects such as 
pigment alteration, blistering, and erythema owing to excessive heat generated in the 
tissue [8,9]. Insufficient delivery of laser light can also be problematic, leading to inef- 
fective outcomes or even paradoxical hypertrichosis [10,11]. To minimize such side ef- 
fects, many studies have been performed to determine the optimal set of laser 
parameters to treat various body sites and skin types of patients [12-15]. 

Even with the right intensity of laser light, however, the actual amount of energy de- 
livered to the target presents local variations when the laser treatment spots are not 
evenly applied [16-18]. The potential threat of side effects from this nonuniformly de- 
livered laser light is serious when we consider that these treatments are often delegated 
to nonphysicians or even nonmedical personnel who lack sufficient training [2,19]. Phy- 
sicians must also practice or train on new laser devices, new applicator tips with differ- 
ent dimensions, and different settings of laser parameters [20]. 

Therefore, it is highly desirable to develop an affordable system to evaluate the opera- 
tor s procedural performance in LHR treatments. In these cases, the spatial patterns of 
delivered laser light energy must be quantitatively accessed, but little attention has been 
focused on this issue thus far. 

The difficulty in implementing such a system is the visualization of the irradiated 
area. Because infrared is used and no prompt marks are left on the treated area after 
LHR treatments, the irradiated spots cannot be easily traced. Recently, a group of re- 
searchers proposed a thermovision camera-based laser- visualizing method [21]. They 
successfully viewed the thermal changes in the skin and quantitatively analyzed the de- 
gree of overlaps and omissions. However, the use of an infrared camera is too expensive 
for general use. Moreover, it is only applicable to the postoperative assessment of per- 
formance accuracy and can be potentially erroneous when any type of skin-cooling 
mechanism is used. Other researchers utilized photosensitive paper and a camera to 
visualize the laser light [22] for investigating intensity profiles across a single spot to 
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analyze the nonideal characteristics of an LHR device. However, this methodology cannot 
be directly applied to the LHR performance evaluation system, primarily because photo- 
sensitive papers are not reusable and extra steps are required to digitize the results. Some 
commercial devices incorporate an auxiliary visible laser light to guide an operator for 
aiming the position of the laser spot; however, these are inapplicable to contact-type LHR 
devices and the trajectory of the irradiated laser spots cannot be traced. 

In this research, a relatively simple system is proposed to visualize and analyze the 
delivery patterns of laser sources during a simulated LHR procedure. The system is 
intended for preclinical uses to evaluate the proficiency of operators and features af- 
fordability and simplicity, based on an off-the-shelf PC camera and digital image pro- 
cessing methods. The purpose of this study is: 1) to demonstrate that the proposed 
system can quantitatively reflect the performance level of LHR treatments, and 2) to 
test the applicability of the system in the field by training novice operators and compar- 
ing two different treatment techniques (sliding versus spot-by-spot). For these pur- 
poses, six novice operators are recruited and their performance is evaluated according 
to four performance indices during the simulated procedures of LHR. 

Methods 

Experimental apparatus 

In the experiments, commercially available laser equipment, a PC camera, and polari- 
zers were used with a simulation bed built in-house. Table 1 summarizes the specifica- 
tions of the devices. 

The simulation bed was made of aluminum profiles in a cuboid shape that was open 
on each side. To mimic the friction of human skin during LHR simulation, a sheet of 
semitransparent silicon rubber was placed on transparent glass on the top of the simu- 
lation bed. Inside the simulation bed, an angled mirror was installed (at 45 degrees) to 
reflect the laser light coming from the top to the front of the simulation bed. Figure 1 
shows the simulation bed with the positioned camera. 

A PC camera (SPC-A30M, Samsung, Seoul, Korea) was placed at a distance of 
200 mm from the front side of the simulation bed to capture the laser light. Its position 
and viewing direction were adjusted to fully view the treatment area reflected in the 
mirror and to not distort the scene. By attaching a series of polarizers (visible linear 



Table 1 Summarized specifications of the devices used in the experiments 



Devices 

Laser 

equipment 

Simulation 
bed 

Silicon layer 
Glass layer 
Mirror 

Camera 



Manufacturer 

LightSheer XC, 
LUMENIS, Inc. 



Polarizers 



Anonymous 

Anonymous 

Anonymous 

SPC-A30M, 
Samsung, Inc. 

Visible linear 
polarizing film, 
Edmund Optics, Inc. 



Size 

20 X 20 mm^ (applicator tip) 
12x12 mm^ (laser window) 

400x250x250 mm^ 

400x250 mm^ 
400x250 mm^ 
400x350 mm^ 

50 X 50 X 70 mm^ (approx.) 



Six orthogonally aligned 
15x15 mm^ film cuts 



Others 

Fluence: 10-100 J/cm^ Repetition: 1-2 Hz 
Wavelength: 800 nm 

Made of 1 5 mm x 1 5 mm aluminum profiles 



To mimic the skin friction 

To mechanically support the silicon layer 

To reflect the laser from the top to the front 

Sensitivity: Visible and infrared regions Frame 
rate: 30 Hz Resolution: 640 x 480 pixels 

transmission: > 40% at 800 nm 
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Figure 1 Simulation bed used in this experiment. Lasers emitted on the human skin layer were 
reflected to an infrared camera. A silicon rubber sheet was stacked on top of the simulation bed to mimic 
the friction of human skin and a normal PC camera was utilized for the visualization of the laser. 



polarizing film, Edmund optics, Harrington, NJ) in front of the camera lens and taking ad- 
vantage of the PC cameras sensitivity to the infrared lights, we could make it function as 
an infrared camera. Additionally, the polarizers protect the camera by significantly attenu- 
ating the intensity of input laser light and increasing the signal-to-noise ratio of the video 
images by removing background noises. Because the intensity of the laser light is much 
higher than visible-band light, the captured video images only contain the laser spot. 
Ideally, two orthogonally aligned polarizers can block all of the incoming lights. However, 
owing to practical discrepancies, six layers of polarizers were used in this experiment. 

Contact-type diode laser equipment (LightSheer XC, LUMENIS, Santa Clara, CA) 
was used in the experiments. The device fires a single pulse of laser light when the trig- 
ger button on the laser applicator is pressed or a train of continuous pulses at a predeter- 
mined rate when the button is held. The fluence and repetition frequency of the laser are 
configurable, but were set to 25 J/cm2 and 2 Hz, the most common settings used in 
clinics, during all experiments. The size of the applicator tip was larger than the size of 
the laser window owing to the cooling area located at its perimeter (Figure 2). 

Study design 

Six novice operators who had no prior knowledge about the LHR procedure were re- 
cruited (all men, 24-31 years) and simulated LHR treatment on the top of the simulation 
bed. Operators were asked to achieve a delivery as uniform as possible of the laser on a 
rectangular target area designated as 140 mm x 90 mm. The usage of gel was mandatory 
to simulate the real conditions of an LHR procedure. Prior to the experiments, a didactic 
lecture was given by an experienced dermatologist about the safety issues, principles, and 
techniques of LHR treatment. A real demonstration followed the lecture. 

Two different techniques of treatments, spot-by-spot (SBS) and sliding, which are 
most frequently used in clinics, were simulated. In the SBS technique, a single laser 
pulse was fired at a time while the laser applicator was repeatedly placed on and re- 
moved from the skin by an operator. In the sliding technique, laser pulses were con- 
tinuously fired while the laser applicator was slid, maintaining contact with skin 
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throughout the whole procedure. The order of the simulated techniques was randomly 
assigned to avoid the possible effect of habituation to the experimental apparatus. 

Each simulated technique of treatment was composed of two separate sessions; the in- 
frared camera installed inside the simulation bed recorded all simulated sessions. To 
simulate the training of novices, feedback on the procedural performance was given to the 
operator between the two sessions. Feedback was provided by viewing the image- 
processed video recorded during the prefeedback session, and erroneous laser delivery 
patterns such as overlapping and omissions of laser spots were confirmed automatically. 
Before initiating each session, operators were allowed free time to test and become accus- 
tomed to their strategy of treatment. Time for the free run was unlimited and varied 
among operators, ranging from 1 to 3 min. Figure 3 shows the overall study design. 

This experimental protocol was approved by a local institutional review board (IRB No. 
C- 1302-075-467 at Seoul National University Hospital) and conducted in accordance with 



SBS technique 




Figure 3 Design of simulation experiments. Two different treatments techniques were simulated and 
eacli treatment teclinique was composed of two recorded sessions of tlie simulation (filled circles). 
Feedback on the procedural performance was given between these two successive sessions. 
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the principles of the 2004 version of the Declaration of Helsinki. All subjects signed the 
informed consent form. 



Image processing 

A laser energy distribution map (LDM) was synthesized through a series of image pro- 
cessing algorithms on recorded video frames. Figure 4 shows the block diagram of the 
overall image processing steps. 

First, calibration between physical and image space was made by locating four cor- 
ners of the target area with 10 laser exposures. The calibration process revealed that 
140 mm of physical dimension was equivalent to 397 pixels in the image space. The de- 
gree of distortion in the camera view and the errors in locating the laser spots were 
inconsequential. 

Next, the frames containing single pulses of the laser were retrieved from the video 
recording of a simulated session. Because a laser spot is shown as a cluster of bright 
pixels in the image, the laser-exposed frames can be selected from a plot of mean gray- 
scale of each image frame. As exemplified in Figure 5, peaks in the plot represent laser- 
exposed frames. Some frames were not counted as valid laser-exposed frames; for 
example, those with no laser spot but an overall increase in brightness level at less than 
30% of normal laser-spotted frames. These frames occurred when the laser was fired in 
the air, mostly during directional changes of the laser applicator during the simulated 
SBS technique. 



Video to frames 



Get laser-exposed 
frames 



In case, more tlian one 
clusters aie found 



Binarization 



Get the cluster with 
max area 



Count the number of 
white pixel-clusters 



Locate the centroid 
coordinates of the laser 
spot 



Calibration info. 



Laser-window template 



Overlay 
the template image 



] 




From a \ideo 
recorded duiing a 
simulation session 



Iteration on each 
laser-exposed fiames 



Reconstructed laser- 
distibution map, 

L(m,n) 

Figure 4 Block diagram of image processing steps to synthiesize an LDM. Laser-exposed frames were 
extracted from the video recorded during a simulation session. Eacli laser-exposed frame was binarized to 
clearly detect the laser spot, then its centroid position was computed. The LDM was reconstructed by 
overlaying the template image of the laser window to each centroid position of the detected laser spot. 
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Figure 5 Mean intensity plot of the recorded video to specify laser-exposed frames. Laser-exposed 
frames were extracted from the video based on the brightness scale of frames. Peal<s in the plot represent 
laser-exposed frames; however, some peaks with less than 30% of normal peak height were not counted 
as valid laser-exposed frames. These frames occurred when the laser was fired in the air, mostly during 
directional changes of the laser applicator when the SBS technique was simulated. 



Next, the position of the laser spot in each laser-exposed frame was located (Figure 6). 
Images were converted to black and white (binarization) with a threshold of 30% of 
maximum value in the frame. Then, the coordinates of the centroid of a white cluster 
were computed. Normally, the converted images contained only a single cluster of 
white pixels. However, in some cases, the reflected laser light from the aluminum pro- 
file appeared as additional white clusters. In such cases, the cluster with the largest area 
was regarded as the true laser spot. 

Finally, a square template that represents the actual shape and size of a laser spot 
was overlaid at the position of the detected laser spot (Figure 6). The gray value of the 
pixels inside the template was set to a constant of 1, which reflects a uniform distribution 



raw images 

locating centroids 
template overaly 

reconstruction of 
laser distribution 
map (LDM) 

Figure 6 Illustration of LDM synthesis by locating laser spots and overlaying template images. The 

raw images of laser spots were binarized and a square template mimicking the actual size and shape of the 
laser-emitting window was overlaid to the centroid position of each laser spot by allowing superimposition. 
The reconstruction process was necessary because distortions were evident in the raw images, caused by 
the scattering effect of the rubber sheet and the point spread characteristics of the camera system. 
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of the energy delivered from the laser light. This type of reconstruction process 
was necessary because the captured images of laser spots were blurred owing to 
the scattering effect of the silicon rubber sheet as well as the point spread charac- 
teristics of the camera system, A calibration process determined the size of the 
template to be 34 x 34 pixels. The LDM was built by accumulating all of these 
template-represented laser spots in the retrieved frames, allowing superimposition at 
overlapped regions (multiple doses of laser). Each pixel in the target area was ini- 
tially set to zero and accumulated its value by the gray value of the template every 
time it was attributed to a laser spot. 

Performance indices 

For the quantitative analysis of the LDM, four performance indices are proposed by 
using general statistics. The first two indices, 5o and 5^, are measures of the errors in 
the LDM, defined as the percent ratio of untreated and redundantly treated areas to 
the total area of the target, respectively, which is shown in Eqns. 1 and 2. 

6o[%]=^^^xlOO (1) 

§z [%] = ' ^ X 100 (2) 
Where 

M N 

Ak = ^^L(m, n)o k k = 1, 2, a and Aq = M x N (3) 

m=ln=l 

f l,if x> y 
1^ 0, otherwise 

(4) 



Here, a is the highest pixel value found in the LDM, L(m,n), which denotes the max- 
imum redundancy in the laser delivery. M and N represent the dimensions of the image 
in pixels (number of rows and columns, respectively). The maximum value of 6^ may 
exceed 100% because the target area can be treated redundantly with more than two 
laser exposures (multiple doses are counted by summing each AiJ. 

The (i is an index that represents an operator s estimation of the spot size, which is 
defined by using the mean of every two consecutive spot distances, dc, as depicted in 
Eqn. 5. 

(i [mm] = mean(dc) x C (5) 
where 

dc= {xi|xi= |S(i)-S(i + l)|,i = l,2,...,(|3-l)} (6) 

The constant C is the conversion ratio between the physical and image spaces, which 
was found to be 0.35 mm/pixel in in the calibration step of our experiments. S is the 
array containing the centroid of the position of each laser spot. |3 is the number of laser 
spots exposed to the target. The ideal value of is equal to the actual size of laser win- 
dow (which is 12 mm in this study). 
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The u is the measure of randomness in the LDM. It is formulated as a normalized 
form of distance variations measured from each laser spot to its nearest one: 

_ Std(dn) 



mean(dn^ 

where 



X 100 (7) 



dn = {Yi I Yi = min( |S(i) - S(j)|), Vje {1,2,...,|3}} (8) 

An increase in u may result in increases of both or either 5o and 5^; however, the re- 
verse is not always true. 

To illustrate the computation of performance indices, a synthetic LDM with three 
shots of lasers is shown in Figure 7. In this example, the size of the LDM and the laser 
spot were set to 13 x 12 and 5x5 pixels and the conversion ratio, C, was set to 1 
(i.e. 1 mm = 1 pixel). The numbers in the pixels indicate the amount of laser ex- 
posure at the site, and the centroid position of each laser spot is marked in red. 
The computation process is as follows: 

1. Maximum redundancy in the laser delivery a = 3. 
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Figure 7 Synthetic LDIW to illustrate the computation process of performance indices. The LDM has 

three overlapped laser spots and the pixel values indicate the amount of laser exposure at the site. Here, 
the maximum redundancy in laser delivery is 3 and the centroid position of each laser spot was marked as 
red. The indices 5o and 6^ are computed based on the pixel values, and \x and u are computed based on 
the centroid position of laser spots. 
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2. Aq, Ai, A2, and A3 are 156, 58, 15, and 2, according to Eqn. 3. 

3. 60 = 62.82% and 6^ = 10.90%, according to Eqns. 1 and 2. 

4. S = {(7,5), (5,8), (9,7)}. 

5. dc = {3.61, 4.12} and dn = {2.83, 3.61, 2.83}, according to Eqns. 6 and 8. 

6. (i = 3.87 mm, u = 0.15%, according to Eqns. 5 and 7. 

Statistical comparisons 

A total of 24 EDMs (six subjects, two techniques for treatments, and pre/postfeedback) 
were collected for the statistical comparison. 

To validate the efficacy of proposed indices in assessing the performance of the LHR, 
an experienced dermatologist sorted the collected EDMs into three groups according 
to performance. Eight EDMs that presented the highest level of performance were se- 
lected first and designated as group G. The other 16 EDMs presented relatively poor 
performance compared to group G, but were sorted again into groups P+ and P_ be- 
cause there were alternative reasons for poor performance. Therefore, EDMs having 
too many overlaps were assigned to group P+ and those with too many omissions to 
group P_, with eight EDMs each. 

The effect of the training was examined by comparing data between prefeedback 
and postfeedback sessions. The comparison was made for each treatment technique 
by comparing Dl and D2 for SBS technique and D3 and D4 for the sliding tech- 
nique. Because we also hypothesized that there is a difference between treatment 
techniques, Dl was separately compared to D3 and D2 to D4, with the presence of 
feedback. 

Statistical analysis of the experimental data was performed by using the Student s t- 
test for the matched paired, one-sample t-test, one-way ANOVA, and a Turkey post- 
test at a significance level of 0.05. Prior to the t-test and ANOVA, the normality of the 
data was assessed by using Shapiro-Wilk s method. 

Results 

Validation of performance indices 

The results of ANOVA indicated that the mean values of performance indices are 
significantly different among performance groups (Figure 8). Group G showed 
lower values of 6^ and u than group P+. The u value of group G was also lower 
than that of group P_. Group P+ and P_ could be distinguished by using any single 
performance index except for u. The one-sample t-test result for (i showed that 
only group P+ presented a significantly different value of from the actual laser 
window size of 12 mm (Table 2), 

To exemplify different levels of performance, nine EDMs selected from three per- 
formance groups are shown in Figure 9. It is clear that the EDMs from group G pre- 
sented better laser distribution than the rest. Specifically, groups P+ and P_ exhibited 
more overlapping and omission, respectively, than group G. In these cases, the mean 
value of 5z was 29.08 for group P+ and 6.05 for group G. The value of was 23.77 for 
group P_ and 13.12 for group G. Therefore, the quantitative analysis results are consist- 
ent with the results of visual inspection. 
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Figure 8 Differences in indices among performance groups. Each performance group was characterized 
by using a combination of indices (T < 0.05, """P < 0.01). The quantitative description of satisfactory patterns in 
laser delivery in LHR treatment aims to achieve minimal randomness in the spatial distribution of laser spots, 
with accurate estimation of the size of the laser spots. 



Training of operators 

Four categories of data described in Table 3 were considered to study the effects of 
training (prefeedback versus postfeedback). The mean values in Table 4 represent the 
differences between Dl and D2 for the sliding technique and between D3 and D4 for 
the SBS technique. 

The paired t-test results of the training indicate that the feedback was effective when 
the sliding technique was used, as indicated by the reduction in 5o and (i. The value of 
(i decreased from 13.49 mm to 11.70 mm, approaching the ideal value of 12 mm. The 
increase in 5^ was not desirable, but the degree of its change was not statistically sig- 
nificant. In contrast to the sliding technique, none of the performance indices changed 
meaningfully when the SBS technique was used. The values in Table 4 represent the 
differences between the indices in the postfeedback session computed relative to the 
prefeedback session. 



Table 2 Deviation of \x from actual laser window size 





H: 12 mm 






(mean ± SD) 


P 


Group G 


+0.50 ±2.1 7 


0.537 


Group P+ 


-1.65 ±2.1 


0.002** 


Group P_ 


1 .59 ± 2.74 


0.144 



'P<0.01 Negative value indicates shorter spacing between consecutive laser spots. 
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Figure 9 Comparison of LDMs thiat represent different levels of performance. Nine LDMs selected 
from three performance groups are shown to exemplify different levels of performance. The LDMs from 
group G presented better laser distribution than the rest. Specifically, group P+ and group P_ exhibited 
more overlapping and omission, respectively, than group G. 



Comparison of treatment techniques 

To compare the two different treatment techniques, the differences between Dl and 
D3 were computed for the case of prefeedback. The same procedure was followed be- 
tween D2 and D4 for the case of postfeedback. The results are summarized in Table 5. 
During prefeedback sessions, operators showed a higher value of u when the sliding 
technique was used than when the SBS technique was used. The remaining indices did 
not show differences between the treatment techniques. The differences were measured 
relative to the SBS technique. 



Discussion 

As pointed out in a previous study [23], methods to learn clinical skills are changing as 
the opportunities for learning through work with actual patients has diminished. The 
same holds true for LHR treatment because current residency programs in dermatology 
place insufficient emphasis on photodermatology and laser therapy [24], and even non- 
physician treatments have become prevalent [25]. The absence of validation methods is 
a problem in LHR treatment; because of this, no operators, even licensed ones, can ini- 
tiate treatment with great confidence. The use of models or simulators has been the 
common practice to tackle this shortage of experience [23] . 

The proposed LHR evaluation tool successfully visualized simulated patterns of laser 
delivery and evaluated them in terms of four indices. Significant differences in indices 
among performance groups were found; however, none of these indices were sufficient 



Table 3 Categories of collected data for statistical comparison 



SBS technique 




Sliding techinique 




Prefeedbacic 


Postfeedbacl< 


Prefeedbacic 


Postfeedback 


Dl 


D2 


D3 


D4 
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Table 4 Effects of feedback depending on treatment technique 




Sliding mode 




SBS mode 






(mean ± SD) 


P 


(mean ± SD) 


P 


6z [%] 


+6.00 ± 8.03 


0.127 


+ 1.07 ±5.27 


0.641 


6o [%] 


-6.32 ±5.1 2 


0.029* 


-4.68 ±6.1 4 


0.121 


M [mm] 


-1.79 ±1.67 


0.047* 


-0.32 ± 1 .08 


0.504 


u [%] 


-0.77 ± 3.83 


0.645 


+0.75 ± 1 .22 


0.188 



Prefeedback |j values were >12 mm, regardless of techniques. 
(13.49 In sliding, and 13.86 In SBS). 

Positive signs Indicate a higher Index value In postfeedback session. 
*P < 0.05. 



to characterize all of the performance groups at once. Therefore, the proposed indices 
should be used in combination to assess the level of performance. 

The ideal situation of group G, the most well-performing group, should have a \i 
value of 12 mm and the lowest values in the rest of the indices. However, only the two 
indices, \i and u, corresponded to this expectation. Therefore, the quantitative descrip- 
tion of satisfactory patterns in laser delivered LHR treatment aims to achieve minimum 
randomness in the spatial distribution of laser spots and accurate estimation of the size 
of laser spots. A moderate level of omissions and overlapping appears to be allowed; 
however, it has been deduced that the degree of overlapping is a more important factor 
than that of omissions when describing satisfactory performance, because group G 
could be distinguished from one of the poorly performing groups in terms of 5^ but 
not Sq. 

The training was effective only when the sliding technique was used, as indicated by 
the improvements in 6o and [i values (Table 4). A careful examination of the changes in 
[I values is given in Figure 10. The more deviation an operator exhibited in [i from its 
ideal value of 12 mm, the more improvements were observed after training. The re- 
duced 6o after the training is believed to be the result of a more accurate [i. In contrast 
to the sliding technique, operators' performance did not significantly improve when 
using the SBS technique. The difference of training effects between treatment tech- 
niques is hypothesized to be largely attributable to the characteristics of the gel. In the 
SBS technique, gels bear the marks of previous contact with the applicator tip at every 
departure of the applicator tip from the target surface. Operators may have been misled 
by the illusion of marks that inhibited them to reflect on the feedback from the system. 
This was confirmed by interviewing the operators as a group after finishing the experi- 
ments. Five of six operators replied that they took advantage of the marks on the gel to 
easily determine the position of the laser spot during the simulation of the SBS 
technique. 

Table 5 Results of LDM analysis on the effects of treatment technique 

Prefeedback Postfeedback 

(mean±SD) P (mean±SD) P 

6z [%] -2.87 ± 8.50 0.445 +2.06 ± 1 2.80 0.71 0 

6o[%] +5.40±6.78 0.108 +3.76± 13.31 0.520 

M [mm] +1.63 ±1.58 0.053 +0.16+1.45 0.794 

u [%] +7.21+5.79 0.028* +5.69 + 5.79 0.146 

Positive signs Indicate higher Index values In the sliding technique. 
*P < 0.05. 
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Figure 10 Individual changes of |j after training for the sliding technique. The greater the deviation 
an operator exhibited in \i from the ideal value of 12 mm, the greater improvement was observed after 
training, except for subject 6. 



In the sliding mode, however, operators were not deceived by marks on the gel, be- 
cause an operator had to adaptively control the amount of force applied to the laser ap- 
plicator tip to cope with the varying frictions while maintaining contact with the 
surface. The results listed in Table 5, which compare treatment techniques, support this 
hypothesis. The higher value of u in the sliding technique indicates that the operators 
experienced more difficulties maintaining constant distances between successive laser 
spots while using this technique, particularly before feedback was provided. This result 
also supports the reason for recommending the SBS technique to novice operators in 
clinics. Therefore, operators are required to be well-informed about the effects of the 
gels on performance, depending on the technique used for the LHR treatment. 

It was interesting to note the trend of elevated prefeedback values, regardless of 
treatment techniques (Table 4). This indicates that the operators overestimated the size 
of the laser spot, which is assumed to be attributable to the mismatch in size between 
the applicator tip and the laser window. Because the actual size of the laser-emitting 
window was smaller than that of the applicator tip, and only the top view of the appli- 
cator tip was visible during the treatment, operators were susceptible to overestimating 
the size of the laser spot. 

The presented study has three limitations, described as follows. First, the aforemen- 
tioned criterion of satisfactory delivery of the laser in LHR treatment can be prema- 
turely generalized, because the performance groups presented in this paper reflect only 
a single professional dermatologists opinion. More general wisdom— for example, a 
composite score of LHR proficiency— will be found when additional tests are conducted 
with a group of dermatologists by comparing their simulated results against those of 
novices. Second, owing to the practical limitations in utilizing clinical settings, the sys- 
tem was tested by a small number of operators; therefore, the statistical analysis 
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performed in this study may have a high probabiUty for type II error. In other words, 
the performance indices might be more significant for characterizing performance 
groups and observing the effects of training. Additionally, a prospective study might be 
conducted to confirm the long-term efficacy of the system in operator training; how- 
ever, the usage of simulators has proven effectiveness in other medical fields for the ac- 
quisition of particularly junior levels of maneuverability [23]. Third, the system might 
not be applicable to noncontact types of LHR devices, so the simulation of treatments 
on curved areas of the body would not be feasible. Because a mirror was used to reflect 
the laser beam to the camera, the location of the laser spot can be erroneous as the 
laser applicator deviates from the vertical line to its target surface. However, this limita- 
tion does not seriously degrade the fidelity of the simulation. Referring to the results 
from a research group that recently analyzed patients' thermal deposit images during 
cosmetic laser procedures [21], the patterns in experts' laser delivery (Figure 7 in the 
study [21]) resembles those of well-performing groups in our experiments (Figure 9 in 
this experiment). In the previous study, operators also tried to cover the target area 
with uniform but not overlapping laser spots, as in our experiments. Although the per- 
cent errors of the omitted and overlapped areas were not directly comparable owing to 
some differences in the methods used for computing indices, the similarities between 
the simulated and clinical trials ensure the applicability of the proposed system in oper- 
ator training. The pulse rate of the laser is another aspect of hardware limitation. A 
normal webcam can capture up to 30 frames of images per second; therefore, lasers 
pulsed higher than this rate may not be correctly detected. However, such settings are 
uncommon in clinical practice. 

Even with these limitations, it is expected that the proposed system could be effect- 
ively utilized in clinics, because it is a cost-effective and intuitive solution to visualize 
and evaluate the proficiency of LHR treatment. The system might be improved to be 
used in other areas of photomedicine. For example, laser applications in the treatment 
of pigmented lesions and facial rejuvenation require different degrees of laser overlaps 
and patterns, depending on the need of the patient and the treatment settings [26,27]. 
Therefore, the performance indices might be adjusted to reflect the general guidelines 
or an expert s strategy for various treatments. An improvement to the software is also 
expected. Accounting for thermal relaxation time, rather than simply counting the 
number of redundant exposures, will be more advantageous in predicting actual ther- 
mal damage to the tissue. 

Conclusions 

A uniform laser delivery during LHR treatment is significant for safe and effective treat- 
ment. A highly affordable system, which is also easy to operate, has been developed in 
this research. The proposed system was able to visualize and evaluate laser patterns 
during preclinical trials without using an expensive infrared camera. For the study, four 
useful performance indices were proposed for assessing the proficiency of operators 
during LHR treatment. With these indices, the developed system could quantitatively 
analyze an operator s proficiency in LHR treatment. A performance analysis of the pro- 
posed, affordable system has shown that operators reduced omission errors by 6.32% 
and accurately estimated the spot distances to match the actual size of the laser-emitting 
window. Further, the proposed system was used as a scientific tool for the comparison of 
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two different treatment techniques (sliding versus spot-by-spot) by observing the perfor- 
mances. Therefore, the proposed training system is expected to benefit many operators in 
clinical practice and to maintain skilled performance in LHR treatment, which may result 
in eventually accomplishing a uniform laser delivery treatment. The proposed system may 
also be applicable to other areas of photomedicine. 
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